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Abstract. Kinetic numerical simulations of the evolution of the Weibel instability during 
the full nonlinear regime are presented. The formation of strong distortions in the electron 
distribution function resulting in formation of strong peaks in it and their influence on the 
resulting electrostatic waves are shown. 



1. Introduction 

The process of magnetic field generation in a plasma is one of the most important problems for 
both laboratory and astrophysical plasmas. Several mechanisms of magnetic field generation 
have been analyzed that are effective in different plasma regimes. At high frequencies, i.e., on 
short electron time scales over which the ions can be assumed to be immobile and collisions are 
unimportant, the Weibel instability can efficiently generate magnetic fields in plasmas with an 
anisotropic electron temperature distribution with zero real frequencies. 

The Weibel instability is an electromagnetic plasma instability driven by the presence of 
temperature or electron momentum anisotropy pQ. 

In this paper we discuss the long time, kinetic, nonlinear evolution of the Weibel instability 
in an electron-ion plasma by using Vlasov-Maxwell simulations in the 1D-2V phase space. 
Relativistic effects on the plasma dynamics will be not considered, which become significant 
only for very high beam speeds. 

The Weibel instability developed due to the counter streaming of electrons have been studied 
and showed that other than the electromagnetic instability, there is also an electrostatic limit 
that corresponds to the well known two-stream instabilites. For more results we refer to 
Refs. [2J[3]. In this paper instead, we assumed temperature anisotropy to study Weibel instability 
and will show, even in this case too, some kind of electrostatic instability develops due to 
acceleration of the particles. 

The Weibel instability emerges from a wide white noise spectrum of wavelengths around the 
electron skin depth, d e . This length plays the role of the natural scale length for the instability 
and will be used in the following as the normalizing length. 

In general, since the growth rate of the Weibel instability is comparable to the electron 
plasma frequency, and also the electromagnetic nature of the instability, this study of the Weibel 
instability in the nonrelativistic limit has been limited so far to the electron dynamics. This 
approximation, which considers the ions as fixed neutralizing background, is valid during the 
linear and very early nonlinear phase of the instability. On the other hand, for longer times, the 
ion response cannot be neglected. 



The paper is organized as follows. In the next section we will introduce the governing 
equations and the initial conditions. The numerical simulations are presented in section III. 
Conclusions follow in section IV. 



2. The Equations 

The kinetic evolution of an electron-ion collisionless plasma is described by the Vlasov equation 
in the phase space, coupled self-consistently to the Maxwell equations. We have used 1D- 
2V (x ,v x ,Vy) phase space, with periodic boundary condition along x direction, unidirectional 
magnetic field B z and a two-component electric field E = (E x ,E y ). 

The Vlasov-Maxwell equations are normalized by the following characteristic quantities: the 
electron mass m e , the characteristic electron density n e , the speed of light c, the electron plasma 
frequency u pe = (4irn e e 2 /m e ) , the electron skin depth d e = c/oj pe and the characteristic fields 
E = B = m e Lo pe c/e. In dimensionless limits the Vlasov-Maxwell equations for electrons read as 

||v||(E|vxB).^ = (1) 
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dB z _dEy dEy _dB ± , 

dt ~ 8x ' dt ~ 8x + y 1 j 

where f e (x,v x ,v y ,t) is the 1D-2V electron distribution function (edf). The Poisson equation 

given by = (rij — n e ), where, </> is the potential and n« and n e are ion and electron density 
respectively. 

We integrate Vlasov equation in the interval L x = 2n/k to have the total length of the box 
large enough to include the most unstable mode in the simulation. The initial dimensionless 
distribution function considered is: 
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fe(x,v x ,v y ,0) = (— Wttt )exp —{-vl-^-vl) + cos(kx + <p(k))] (3) 
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where, L x is the length in the x direction, v t ^ is the thermal velocity of the electrons and 
T x and T y are the temperatures along x and y direction respectively, the modes of the wave, 
k = 2irn/L x , and e is the perturbation amplitude. To include all the modes possible we choose 
300 as the summation limit. At t = we introduce also a perturbation on the magnetic field: 

300 

B z = a in ^2 [cos(k n x + (p(k))} (4) 

n=l 

where, a« n is the magnetic field amplitude and is independent of the mode. 
3. Numerical Simulations 

We have performed the numerical simulation (see Ref. [1] for the numerical algorithm) of the 
Weibel instability in plasma, by integrating numerically Eqs.(l)-(2) with the initial conditions 
given by Eqs.(3)-(4). In particular, the electron and proton distribution function is advanced 
in time using the electromagnetic splitting method in the 1D-2V limit. The Maxwell equation 
are integrated in time by using the third order Adams - Bashforth algorithm and standard Fast 




Figure 1. Fig(A) shows the Fourier transform of the magnetic field vs k at t = 225 showing the 
most unstable mode at k = 1.25. Fig(B) plots the time evolution of .B^fc (solid lines), £ Xi fc(doted lines), 
.E^k (dashed lines) for k = 1.25. 



Fourier Transform technique are applied to calculate the space derivatives. The temperature 
anisotropy considered is T y /T x = 12. We have performed the simulation with higher and 
lower temperature anisotropies also. But with lower temperature anisotropies the growth of 
the magnetic field is not strong enough to see the features we are interested in and also take 
longer time to run the simulation to reach the nonlinear regime. Nothing interesting is found at 
higher temperature anisotropies too. Hence, for the sake of qualitative agreement we assumed 
the particular case of anisotropy with T y /T x = 12. The thermal velocity is taken as vth = 0.02, 
corresponding a plasma with a non-relativistic temperature. The perturbation amplitude for 
each wave is ai n = 1(T 4 . Initially, we consider a field free plasma, hence < B z >= 0. The 
number of points in the x and v directions are N x = 600 and N Vx = 90, N Vy = 90. The time 
step is At = 0.005. The evolution of the system is investigated up to t = 480w~ e 1 . The maximal 
velocity which defines the velocity interval is set to v Xmax =0.1 and v ymax = 0.3 in the x and y 
directions respectively such that —0.1 < v x < 0.1 and —0.3 < v y < 0.3. The length of the space 
domain is L x = 6ir. The space discretization in x is dx = 0.03d e , where d e is the electron skin 
depth. 

In Fig.l (A) we plot the Fourier amplitudes of the magnetic field B z ^ vs k at t = 225. We see 
the modes in the interval 0.1 < k < 3 have grown with a maximum growth rate corresponding 
to k = 1.25. Note that a small portion of interval in k is taken. In (B) we plot the Fourier 
amplitude of the magnetic field and of the x and y components of the electric field. In the 
linear phase, t > 100, the magnetic field grows exponentially, corresponding to a growth rate 
7 = 0.0390. During this phase the Fourier amplitudes of the electrostatic field E x ^ remains 
negligible. However from t = 190 to t < 300, the electrostatic field starts to grow at a very 
rapid rate, of the order of twice the linear growth rate of the instability, 7 es = 27(7^ is the 
electrostatic growth rate), a clear signature of the beginning of the non linear regime. At later 
time, t > 300 the magnetic field saturates. The electrostatic field levels off at t > 300 due to 
particle trapping in phase space. 

In Fig.2 we plot the electron distribution function(edf) in (v x ,v y ) space at time t = 
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Figure 2. Contour plot of the distribution function at t = 0, 300, 320, 380 at x = 7.8, 6.4, 10.2 
which are the points of maximum, zero and minimum amplitude of the magnetic field. 



0, 300, 320, 380. Frames in the first row show the edf at t = at the position where the magnetic 
field is maximum, zero and minimum. Next rows show the same at different times. At t = 
we can see that the velocity in the y direction of the particles is greater than in the x direction, 
corresponding to our initial condition where, in the y direction the particles are hotter in the x 
direction. A deformation of the edf is observed from t > 300. We see that there is a x-dependent 
rotation of the edf at the positions of maximum and minimum magnetic field at each time. On 
the contrary there is net acceleration of electrons in the direction of the electrostatic force at the 
position of zero magnetic field. We also observe that with increase in time additional structures 
in velocity space are formed elongated along v x . These additional structures in phase space may 
correspond to the particles that have been accelerated by the induced electric field and deflected 
by the magnetic field. To see further, we plot in Fig. 3 in frames (A) and (B), the induced electric 
field vs time at x = ir/2k and x = 3ir/2k, positions corresponding to first and second zero's in 
the magnetic field. We see that the induced electric field at first zero is positive while for the 
next it is negative, which relates to the acceleration of the particles as mentioned earlier. Same 
kind of distortions are observed in the corresponding contour plots in (C) and (D). 

As shown earlier, for t < 300 the electrostatic field grows exponentially. In Fig. 4, frames (A)- 
(B)-(C) show the time evolution of the electron distribution function versus v x at fixed spatial 
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Figure 3. Frames (A, B) plot induced electric field E y vs time at x = 6.4, 9.4 where B z = 0. 
Frames (C, D) plot contour of edf in (v x , v y ) phase plane at the same positions. 



(A) t = 




(B) t = 300 




-0.10 -0.05 o.oo o.os 0.10 



-0.10 -o.os o.oo o.os 0.10 



(C) t = 370 




-o.ro -o.os o.oo o.os 0.10 



(D) k=25 




400 



Figure 4. Frames (A,B,C) show the time evolution of edf versus v x at t — 0, 300, 370. Frame (D) plots 
the time evolution of E x ^ at k = 25. 



position. We observe bumps along v x direction of the distribution function at t > 300. These 
bumps remain for quite a long time up to t = 370 with a resonant velocity v r = 0.04. In the 
last frame we show for corresponding k=25 the electrostatic field in time. We observe in the 
plot that indeed the electrostatic field grows in time during the period when we observe bumps 
in the edf. 



4. Conclusions 

With the help of the numerical simulations, we have shown the kinetic saturation of the Weibel 
instability. Furthermore, it has been shown that the evolution of the Weibel instability leads to 
distortions in the electron distribution function that give rise to the resonant onset of electrostatic 
field. 
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